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I. Introduction 



The direct analytical or numerical calculations of the higher-order terms to the physical 
quantities in the concrete renormalization schemes provide the important information about 
the behaviour of the corresponding perturbative approximations. However, there are also some 
other approaches to treat the problem of the extraction of certain information from the trun- 
cated perturbative series. These approaches are the principle of minimal sensitivity (PMS) ^ 
and the effective charges (ECH) prescription [Q], which is equivalent a posteriori to the scheme- 
invariant perturbation theory 0. Of course, it is better to use these approaches directly in the 
concrete orders of the perturbation theory, as was done in QCD in Refs. However, if one 

adopts the point of view that these methods really pretend to the role of "optimal" procedures 
in the sense that they might provide better convergence of the corresponding approximations in 
the non-asymptotic regime, it is possible to try to go one step further and apply the procedure 
of re-expansion of the "optimized" expressions in the coupling constant of an initial scheme. 
One can consider the residual (A^+l)-th order term as the estimate of the (A^+l)-th order 
correction in the initial scheme 

The re-expansion procedure was already applied for the analysis of the perturbative predic- 
tions for — in QED |^ (for related considerations see Ref. and for the estimates of 
the QCD corrections to definite physical quantities. In these works, the quantities under study 
are the Drell-Yan cross-section at the 0(a^)-level |TT|, R{s) = <7tot{€.~^€.~ — > hadrons) / a{e^e~ 
fi~^fi~), Rr = r(r — >• z/^ + hadrons) /V{t ^^17^6), non-polarized and polarized Bjorken sum 
rules at the 0{a'^) and even 0(af)- levels [0, |T3[ and the singlet contribution to the Ellis- Jaffe 



sum rule at the 0(ag) -order 

It is clear that the re-expansion formalism, which is similar to the procedure used in Ref. 
151] to predict the RG-controlable /n(m^/me)-terms from the expression for {g — 2)^ through 



the effective coupling constant a(m^/me), correctly reproduces the RG-controlable terms [jT| 



1^]. One can also hope that it can give the impression about the possible values of the 
constant terms as well. This hope is based on the observation of the existence of a satisfactory 
agreement of the results of application of the re-expansion procedure in QED M and QCD 



n, O] with the results of the explicit calculations. It should be stressed that on the contrary 



to the RG considerations of Ref. |jT5[, the "optimization methods" are dealing with the full 



RG-invariance of the quantities under consideration, which produce the additional equations, 
relevant to the freedom in the choice of the higher-order coefficients of the /3-function. The 
solution of these equations gives the possibility to define the sets of scheme-invariants which 
are the cornerstones of the "optimization" methods. 

However, in the definite cases the procedure of re-expansion of the "optimized" results 
can run against some reef, which was overlooked in the process of some previous applications 
[|I], ^ In the case of the analysis of the perturbative series for {g — 2)^ this problem grows 



out from the non-careful treatment of the light-by-light scattering graphs with the electron loop 
coupled to the external photon line. 

In Sec. I of this work we repeat the decription of the basis of the formalism used by us. The 
exact expressions for the terms in the re-expansion formulas are derived. It is demonstrated 
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that the estimates obtained using the re-expansion of the ECH expressions are identical to the 
results of calculations of the (A^+l)-th order corrections in the special scheme, where all lower 
order coefficients of the physical quantities and the QCD /3-function are defined in a certain fixed 
scheme (in the case of QED the on-shell (OS) scheme is usually used) and the (iV-l-l)-th order 
coefficient of the /3-function coincides with the (A^-l-l)-th order scheme-invariant coefficient of 
the ECH /?-f unction /3e//. 

In Sec. II, using the information about the four- loop coefficient of the QED /3-function in 
the OS scheme Jl^ we generalize the considerations of Refs. |^, |10] to the five-loop level. We 



follow the proposals of Ref. and consider the light-by-light scattering graphs mentioned 
above separately in our RG-inspired analysis. We show that this empirical improvement leads 
to more satisfactory and thus more reliable estimates of the five-loop contributions to {g — 2)^ 
than in the case of non-separation of the light-by-light scattering contributions. Finally in the 
Appendix we present the expressions for the six-loop RG-controllable contributions to the muon 
anomalous magnetic moment which follow from the analysis of Sec. II. 

II. The Description of the Formalism 

Consider first the order 0{a^) approximation of a Euclidean renormalization group invariant 
quantity 

N-l 

Dn = doa{l + J2 dia') (1) 
1=1 

with a = as/iT being the solution of the corresponding renormalization group equation for the 
/3-function which is defined as 

da 

/z^— = Pia) = -Poa'il + E c.aO . (2) 

The coefficients di,i > 1 and Ci,i > 2 are scheme-dependent. In order to calculate them in 
practice it is necessary to specify the scheme of subtractions of the ultraviolet divergences. In 
QED the OS scheme is commonly used. However, this scheme is not the unique prescription for 
fixing the RS ambiguities, which affect the values of these coefficients. In both phenomenological 
and theoretical studies other methods are also widely applied. 

The PMS fl]] and ECH prescriptions stand out from various methods of treating scheme- 
dependence ambiguities. Indeed, they are based on the conceptions of the scheme-invariant 
quantities, which are defined as the combinations of the scheme-dependent coefficients in Eqs. 
(1) and (2). Both these methods pretend to the role of "optimal" prescriptions, in the sense 
that they might provide better convergence of the corresponding approximations in the non- 
asymptotic regime, and thus allow an estimation of the uncertainties of the perturbative series in 
the definite order of perturbation theory. Therefore, applying these "optimal" methods one can 
try to estimate the effects of the order 0(a^^^ )-corrections starting from the approximations 
D°^^{aopt) calculated in a certain "optimal" approach [ll], 0, [T^ . 

Let us follow the considerations of Ref. re-expand D'^^[aopt) in terms of the coupling 
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constant a of the particular scheme 



DTM = DM + dDfa''^^ (3) 

where 

5Dt = q) - , cf ) (4) 

are the numbers which simulate the coefficients of the order 0(a^'^^)-corrections to the physical 
quantity, calculated in the particular initial scheme. The coefficients Qn can be obtained from 
the following system of equations: 

or 

^{Dn + fiiva^+') = 0(a^+2)^ ^ > 2 (5) 

OCi 

where the parameter r = /5o^n(/i^/A^) represents freedom in the choice of the renormalization 
point IX. The conventional scale parameter A will not explicitly appear in all our final formulas. 
The system of these equations can be solved following the lines of ref . [|[] . Let us stress again 
that the difference between the "optimization" equations and the RG-approach of Ref. [|l^ lies 
in the fact that the latter one is dealing with the ffist equation from the the system of Eq.(^ 
only. The quantities Vti can be related to the scheme invariants pi in the following way: 

pi = di + Y^ci - ^li{di, di-i] ci, (6) 

Note that the general expressions of the scheme-invariants pi and of the correction terms Qi can 
be defined in different ways. Various definitions differ by scheme-independent constant terms. 
We are choosing these correlated constant terms imposing the condition that the expressions 
for the scheme-invariants pi are connected with the coefficients cf"^ of the ECH /5-function 

PeffiaECH) = -poalcH ( 1 + ClttECH + (^f^^^^Ecn) (7) 

^ i>2 ' 

as 

^ECH 



pi 



l-l 

where 

D{aECH) = doaEcnia) . (9) 

The concrete expressions for the invariants pi and thus for the correction terms Qi can be 
derived from the following equation: 

PeffiaECH) = (10) 

We present here the final expressions, which are already known |jl|: 

1^2 = dodi{ci + di), (11) 
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= dodi{c2 - ^cidi - 2dl + 3d2) . (12) 
and the new terms which we evaluated 

= ^(Scgrfi + 02^2 - 4c2rf? + 201^1^2 " cic/3 + 14^1 - 28d?rf2 + 5dl + 12^16/3) (13) 



^5 = — (4c4(ii — Scsd^ + 2csd2 — Ac2did2 + 8c2dl — 2cid4 + Qcidi — 16cidld2 + 801^2 

+8cidid3 - 48dl + 120dld2 - 48didl + 16^2^ - SGrf^c/s + 20^1^4). (14) 

These terms reproduce the RG controllable logarithmic contributions. In the case of the five- 



loop level one can reobtain the QED results presented in Ref. |19|. We discuss this point in 
more detail in the next Section. 

It should be stressed that in the ECH approach df^^ = for alH > 1. Therefore one gets 
the following expressions for the higher-order corrections in Eq. (3): 

dD^^'"" = n2{di,ci) (15) 

dDi^"" = ns{di,d2,Ci,C2) (16) 
5Df^^ = fi4(c?l,rf2,C?3,Cl,C2,C3) (17) 

where Q2, ^3 and are defined in Eqs. (11), (12) and (13) respectively. 

One can understand from Eqs. (6), (8) that the expressions for Qn and for the corrections 
5D^^^ in Eqs. (|15D-(17) are the exact numbers which are related to the special scheme. This 
scheme is identical to the initial scheme at the lower order levels and is defined by the condition 
cn = cff^ at the (A^-l-l)-order, where c%^^ is considered as an unknown number. This means 
that the correction coefficients dD^ are related to the initial scheme only partly. However, it 



was shown in Refs. ||T2|, [T3| that in certain cases the numerical values of these coefficients are in 
satisfactory agreement with the results of the explicit calculations. A posteriori we consider this 
fact as an argument in favour of the possibility of the application of the re-expansion procedure 
in the cases discussed by us. 

In order to find similar corrections to Eq. (3) in the N-th order of perturbation theory 
starting from the PMS approach |]l|, it is necessary to use the relations obtained in Ref. pO[ 



between the coefficients rff*^"^ and cf^^ (i > 1) in the expression for the order 0{apj^^g) 
approximation D^^^ {apMs) of the physical quantity under consideration: 

d^""' = + OiapMs) (18) 



i + l\ N 

where c^^^^ = Ci. Using now Eq. (|T8|) one can find the corresponding coefficients of the NLO 



approximation 02^^^^ {apMs) 

d^i""' =-\^^ + 0{apMs) (19) 
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and the related expression for the NNLO correction dDJ^^^^^ 

5D^^' = 5D^^'' + ^ (20) 

where 50^^^ is defined by Eqs. (15), (11). 

Repeating now the similar considerations at the NNLO level we get from Eq. (0) the 
following expressions for the NLO and NNLO coefficients df and (if 

dfA^^ = o + 0(apMs); 

dl''' = -lc^''' + 0{apMs). (21) 

Substituting now Eqs. (pT]) into Eq. (12) one can observe that the corresponding next-to- 
next-to-next-to-leading order (N'^LO) correction 6D^^^^ in the re-expansion formula of Eq. (3) 
identically coincide with 5D^^^ defined by Eqs. (16), (12), namely that 

5^3^^^ = 5^3^^^. (22) 

A similar observation was made in Ref. [^] using different (but related) considerations. In 
fact this expression means that the order 0{apMs) correction to d^^^ is cancelling the leading 
order term in the expression for (if We have checked this feature explicitely. 

In fourth order of the perturbation theory the additional contribution to 6D^^^^ has a more 
complicated structure. In order to get it, it is necessary to substitute the following expressions 

= U + 0{a) 



= -^cr^ + 0(a) 

= -^cr^ + 0(a) (23) 

into the expressions for the scheme-invariants p2 and ps and then into the analytical expression 
for The expression for fl4^{df'^^ ,cf'^^) in Eq. (4), which results from this analysis, reads: 

(in r 1 PflfQ 4 , pfl^Q^Q 5 2 PMS ^ 41 



n4(c^r ^ cm = j[-^c,cr' - ^(cr^)^ - ^c^cf- + —ct] (24) 



where 



2 gv 2 3g i; 



9 7 

-{d2 + C2-dl-cidi + —cl) + 0{a) (25) 
8 35 



and 



cf*'^^ = 4(4 + Jcg - csdi - 3did2 + 2d\) + \ci{d2 + C2 + 3^? - cirfi + -^cl) + 0(a). (26) 
2 2 lUo 



The expressions for Eqs. ( |2^ - (^) are the pure numbers, which do not depend on the choice 
of the initial scheme. We will show in the next Section that in the case of the consideration 
of perturbative series for {g — 2)^ the numerical values of n4((if cf *^'^) are small and thus 
the a posteriori approximate equivalence of the ECH and PMS approaches is preserved for the 
quantities under consideration at this level also. 

In certain considerations we will need to use a generalization of the expression for Q2 to the 
case when the intitial perturbative series is starting from corrections of order 0{a^) with p > 1 



D^P =doaPil + J2d,a''). (27) 

i>l 



In this case the expression for the corrections terms read 

{n^i''^)ECH = ^dodl + c^orfici. (28) 
2p 

The corresponding correction related to the PMS-improved expression was originally obtained 
inRef. |T]]: 

{n^i^)pMS = ^^dodl + dodici + —^—docl (29) 
2p 2{p+lj 



III. Applications to (g — 2] 



It is well-known that the expressions for anomalous magnetic moments of the electron 
ae = {g — 2)e/2 and muon = {g — 2)^/2 are known at the four-loop order from the results of 



calculations of Ref. |2T|] and Refs. [^], respectively. The three-loop correction to is now 



known with more accuracy than previously . Combining the currently available information 
about the coefficients of the perturbative series for Oe and we have the following expressions: 

Oe = 0.5a - 0.3294789.. .a^ + 1.17619(21)a^ - 1.434(138)a^ (30) 

a^-ae = 1.09433583(7)0^ + 22.869265(4)a^ + 127.55(41)a^ (31) 

where the expansion parameter a = a/ir is related to the fine structure constant a and the last 
term in Eq. ( |3lD is the result of the most recent calculations of Ref. |^ stimulated by the 



work of Ref. [Tj]. Combining Eq. (pO]) with Eq. (|3TD we arrive at the following approximate 
expression for a^: 

= 0.5a + 0.76585a2 + 240^ + 126a^ + 0{a^). (32) 



The order 0{a^) correction to is only partly known [^]. Our aim will be to try to touch 
the existing uncertainty due to the totally non-calculated order 0(a^)-contribution to Eq. (|32D 



using the re-expansion procedure outlined in the previous section , which is compatible with 
the RG-formalism. 

It is known that in the OS-scheme the coefficients of the corresponding perturbative series 
depend on the large ln(m^/me)-contributions starting from the two-loop level. The parts of 



these effects are governed by the RG-method E5| (for a recent application of the RG method 
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to a^, see Refs. [T^). However, there are also certain ln(m^/me)-contributions, which 
are not governed by the RG-method. They are associated with the hght-by-hght-scattering 
electron loop insertions coupled to the external photon line. These contributions appear first 
in the three-loop graphs, which were subsequently calculated numerically in the works of Ref. 



2^, ^ and recently evaluated analytically in the work of Ref. 



In view of the different origin of the lower ln(m^/me)-contributions we divide all diagrams 
into two classes. The first class contains all diagrams with an external muon vertex and dressed 



interenal photon lines (see Fig. 1). As well as in Ref. we will not include the diagrams 
with electron loops to which four internal photon lines are attached. However, we will include 
four-loop diagrams typical to but with substitution of the external electron vertex to the 
muon one. The second class of the diagrams includes the diagrams with electron light-by-light 
scattering subgraph, to which three and four internal photon lines are attached (see Fig. 2). 
Let us stress that all log(m^/me)-terms of the diagrams contributing to the first class are 
totally controlled by the RG-method, while in class (II) only the parts of these contributions 
are governed by the RG-technique. 

In accordance with our classification we represent the expression for in the following form 

a, = a^P + a^i'\ (33) 
The concrete contributions to Eq.(33) read 

a« = d^i^ail + d^pa + SPa^ + Spa^ + ...) (34) 
af) = 4^V(l + 4'')a + df)a2 + ...). (35) 

Note, that the coefficients di {i > 1) contain the RG-controllable ln(x) = ln(m^/me)-terms. 
Indeed, the corresponding contributions to are governed by the RG-equation 

where j3{a) is the QED /3-function in the OS-scheme, which is defined as 

da 



m 



' /3(a) = 5: (37) 



To our point of view, the separation of all diagrams to the two classes mentioned above is 
respected by the property of the RG-invariance. At least we do not know the arguments why 
the sum of the diagrams which belong to the class (I) and to the class (II) should not obey the 
RG-equations seperately. 



The coefficients of the /3-function are known at the four- loop level |T^. They have the 
following form 
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/33 



288 
5561 

5184 ~ 
-0.571. 



23 



gC(2) + -C(2)ln(2)--C(3) 



(3^ 



Thus, the related coefficients q = A//5o {i > 1) read ci = 3/4, C2 = —1.26, C3 = —1.713. Let 
us write down the asymptotic expansions of the coefficients of the contributions as 



d 



B2 + C2 Hx) 

Bi + d ln{x) + Di ln^{x) + \n^{x) 

B5 + C5 ln(a;) + ln'(a;) + \n%x) + F5 ln^(a;) 



(39) 



and 



"1 



5i 

52 + C2ln(x) 

53 + C3ln(x) 



(40) 



The coefficients Ci, Di, Ei, Fi and Cj, L'j can be related to the coefficients of the /3-function using 
either the RG-considerations of Refs. W^, ^^, ITyi or the explicit expressions for the coefficients 
Qi and Q^'''^ in the corresponding re-expansion formulas (see Eqs.(ll)-(14) and Eq.(28)). The 
results of the corresponding analysis have the following form 

C2 = 2(3oBi (41) 



C4 
E4 



C; = 4(3oB2 + 2(3,Bi 

D, = APlB^ 

= 6PoB3 + ^PiB2 + 2P2Bi 

= 12/3252 + 10/3o/?i5i 

= 8P|B^ 



(42) 



(43) 



8P0B4 + 6P1B3 + AP2B2 + 2P3B, 
2A(3lB3 + 28/3o/3i52 + QI3lB^ + 12I3^(32B^ 
104 

?>2(3lB2 + —(3l(3,B, 
16/3o'5i 



(44) 
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C2 = ePoBi 

= 24/?2Si. (45) 

Note, that in the case of the diagrams of set (II) the corresponding coefficients Bi, B2 and B^ 
contain the contributions of the non-controllable by the RG method ln(x)-terms. 

Let us first discuss the applications of the procedure of Sec. II to the diagrams of set (I). 
In this case the correction terms Q2 — ^4 reproduce all ln(a;)-contributions presented in Eqs. 
(39). Moreover, one can get from re-expansion procedure the exact values of the constant 
terms Bi {i > 3) which do not depend on the ln(x)-terms. In the case of the application of the 
ECH-improved variant of the OS-scheme these constant terms are defined by the conditions 

D _ O {^OS jOS jOS „ „OS\ 
J->i — i'i-ll,"o '1 ) •••) "4-25 ^^l' "^i— 2/ 

= fi,_i(5f^...,i?,_i,ci,...,cf4). (46) 

Similar terms which arise from the PMS-improved expressions can be obtained after taking into 
account the additional scheme-independent contributions derived in Sec. II. We will demonstrate 
that the numerical values of these contributions in the cases considered by us are not large. 

The concrete values of the coefficients Bi, B2'^ , B^^ are known from a comparison of the 
results of the RG-inspired analysis with the results of the analytical and numerical calculations 
P^]. The coefficient Bi = 0.5 is of course well known. The asymptotic expression of the 
coefficient B2, derived in the limit mg/m^ 0, can be found in Ref. ||2^: -B^'^ = — i| + a[^^ 



-1.022923. The value of the coefficient B^^ = 2.1 Al was obtained in Ref. 12^ after subtracting 
the contributions of the light-by-light scattering graphs of the set (II) and of the RG-controllable 
contribution of Eq. (39) from the expression for the three-loop correction to a^. 

The value of the coefficient B^^ , which will be used by us, is different from the one given 
in Ref. p^. The difference comes from the fact that contrary to the classification of Ref. 



2J], we are including into the considered set (I) the four- loop diagrams typical to but with 
substitution of the electron vertex and internal electron loops to the muon ones. Moreover, it 
is necessary to modify the value of B'^^ presented in Ref. p2| in accordance with the results of 
the analytical and numerical |^ re-calculations of the diagrams with three-loop insertion 
into the internal photon line of the lowest order contribution to a^. 

In order to determine the value of the coefficient B^^ in our case we used the following 
expression 

^ = - Af{^^) - C4 Hx) - D4 In'(x) - hi\x) (47) 

where C4, Z^4 and are determined by Eqs. (43) and the value of ^^^^(77) f» —116.7 is the 
sum of the eight-order contributions of the diagrams with electron light-by-light scattering 
subgraphs p2|. The numerical value of the coefficient B^^ is thus B*^^ = —7.74. 



In order to study the predictive abihties of the re-expansion procedure described in Sec. II 
we present in Table 1 the numerical results of our estimates of the coefficients Bi [i > 3) and 
compare them with the exact results for Bf^ and Bf^ presented above. 
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Order 


BO, 


B^{ECH) 


B^{PMS) 


i=l 


0.5 






i=2 


- 1.022923 






i=3 


2.741 


1.326 


1.396 


i-4 


-7.74 


-5.48 


-5.48 


i=5 




41.6 


41.7 



Table 1: Estimated values of the coefficients Bi for the diagrams of set (I). 

One can see that the re-expansion procedure used by us is reproducing well enough the values 
of the coefficients B^ and i?4 (it gives the correct sign and predicts the order of magnitude of 
these coefficients). Therefore, we hope that the estimate of the five- loop constant term B^ is 
also rather realistic. Notice also the sign-alternating character of the results of the estimates 
presented in Table 1. This feature has something in common with the expectation that the 
RG-improved QED series for the Euclidean physical quantities should have sign-alternating 
behaviour |^ . 



Taking now into account the numerical value of the RG-controllable terms in Eqs.(39), (44) 
we arrive at the following estimate of the five-loop contributions of the diagrams of set (I) into 
a,. 



,(10) 



{I) 



B^^ + 8.55 
bQ.l{ECH) 
50.2{PMS). 



(45 



This estimate is almost non-sensitive to the concrete realization of the method of optimization. 
Notice also the effect of reduction of the value of the RG-controllable five-loop contributions 
presented in Refs. [|^, [l^. Let us stress again that this fact is explained by necessity of the 

r the constant term Bf^ derived in Ref. p2 
First 



modifications of the results used in Refs. 



These modifications come from two ingredients. First, it is necessary to use the corrected 
expressions obtained in Refs. [O, EBI of certain four-loop graphs contributing to and second 
to add to the values of B^^ cited in Ref. p2[ the constant terms due to the four-loop graphs 
typical of but with a substitution of the electron vertex and internal electon loops to the 
muon ones. As is known from the results of Ref. |^ the addition to the considerations of 
Ref. [IS] of the diagrams with the internal muon loops leads to strong cancellations. The 



comparison of the RG-controllable expressions of Eq.(48) with the similar one derived in Refs. 



1^, indicates the same pattern. 



Let us now discuss the applications of the outlined procedure for the estimates of the five- 
loop contributions of the diagrams with the light-by-light scattering subgraphs of Fig. (2). The 
special feature of the application of the re-expansion procedure to the diagrams of set (II) is that 
the corresponding terms B^ depend from the ln(a;)-terms non-controllable by the RG method. 
The most precise value of the coefficient c/q^^^ = Bi = 20.94792... is known from the results of 
the analytical calculations of Ref. [EH] . The numerical result for the sum of the corresponding 



four-loop graphs reads |22 



116.7. 



(49) 
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Using now Eqs. (28), (29) we arrive at the following numerical estimate of the sum of the 
corresponding five-loop graphs 

a(^io)(jj) _ ^(3) _ 520.8{ECH) 

= 525.2{PMS) (50) 

which includes the contribution of both RG-controllable and RG-non-controllable ln(x)-terms. 

The estimates of Eqs. (50), (48) should be compared with the one given in Ref. 

aj^io) = 570(140) (51) 

where the central value comes from the exact calculation of the contributions of the diagrams 
of the set of the light-by-light-type diagrams with two one-loop electron loops inserted into the 
internal photon lines (see Fig. 3 ) and the error bar ±140 stands for the estimate of other con- 
tributions (mainly RG-controllable ones). Our estimate of Eq. (50) is in very good agreement 
with the central value of the estimate of Eq. (^l]), while the estimate of Eq.(48) lies within the 
range of the careful estimate ±140 of other contributions. 

Note, however, that our total estimate of the considered tenth-order terms 

aj^io) = 570.9{ECH) 

= 575A{PMS). (52) 

also includes the contribution of the tenth-order diagrams depicted in Fig. 4 and not included in 
the estimate of Eq. ([511). These diagrams are formed by the insertion of the two- loop electron 
loop into the internal photon line of the lower light-by-light-type diagram. The contribution 
of this set of tenth-order diagrams was estimated in Ref. |]^. In order to understand the 



uncertainties of this estimate better it is useful to write down a RG-relation analogous to Eqs. 
(43) for this set of diagrams separately. Notice, that this contribution should be proportional 
to the two-loop coefficient of the /5-function (which is determined by the graphs inserted into 
the internal photon line). Using this observation we arrive at the following relation 

a'il'^\FigA) = B^{FigA) + QBi(3i \n{x). (53) 

The main contribution to the estimate of Ref. comes from the ln(x)-term. Indeed, it has 
the following numerical value QBi(3i ln(x) = 167.47. This expression should be compared with 
the estimate a'-^^^FigA) = 176 ± 35 given in Ref. [pill. One can see that this estimate is 



relevant to the RG-controllable contribution only. However, from the re-expansion procedure 
we can see that the contributions non-controllable by the RG-methods might be non-negligible 
(see Eg. (0)) and might affect the final numerical value of the diagrams belonging to this set. 
In order to study this guess in detail it is of interest to calculate the diagrams of Fig. 4 explicitly. 
This calculational project is rather realistic 



It is also interesting to understand deeper the uncertainties due to other diagrams which are 
included neither in the "optimized" estimates of Eqs. (48), (50) nor in the original estimates 



of Ref. |22|]- These diagrams, depicted in Fig. 5, form a new class of diagrams, which cannot 
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be touched by the RG-inspired analysis. Indeed, one can hardly expect that any resummation 
procedures dealing with light-by-light-type graphs with three internal photon lines will be 
able to give the estimate of the light-by-light-type graph with five internal photon lines. The 
expressions for the ln(a;)-terms the non-controllable by the RG-method for this type of graphs 



can be read from the considerations of Refs. [^]. The result was used in Ref. where the 
following estimate of the diagrams of Fig. 5 was presented 



a^^^\Fig A) = 185 ±85. (54) 



Combining our estimates of Eq.(0) with the ones of Eq.(|5^) we get the final result of appli- 
cations of the re-expansion procedure supplemented by the estimates of the diagrams of new 
structure which are non-touched by this method 

a^^°) ^ 700. (55) 

Let us stress again that the new ingredient of our analysis, which distinguishes it from 
previous applications of the re-expansion procedures in QED |^, |10|, is the separation of the 
considered initial diagrams to two classes, one of which consists of the diagrams, relevant to 



the effects of "new physics" , discussed in more detail in Refs. [p4| , |33[| . This procedure finds its 
support in the theoretical considerations of Ref. |]TB . 



Moreover, we checked that in spite of the good agreement of the application of the re- 
expansion procedure to the non-separated sixth-order expressions for with results of the 
eighth-order calculations 0, the straightforward application of Eq.(13) to the non-separated 
eighth-order approximation of Eq. (32) results in the non-confortably large tenth-order estimate 
a^^'^^ ~ 2160. It is possible to understand that the reason of the success of the application of 
the re-expansion procedure to the non-separated sixth-order approximation is connected with 
the fact that the use of Eq.(12) (and more definitely its last term) gives for the eight-order 
light-by-light-type term the following estimate a|f^ (77) = 8^0^1^2(77) = cl^uHii) which is 



known to be in good agreement with the results of direct numerical calculations . However, 
at the next level of perturbation theory the expression for the correction term of Eq.(13) 
has a more complicated structure and thus the resulting non-separated estimates turn out to 
be non-confortably large. Moreover, we consider the satisfactory agreement of the results of 
separated estimates with the estimates given in Ref. as an argument in favour of treating 



the diagrams with the electron-loop light-by-light scattering graphs separately. 

Another interesting question is connected with the problem of the comparison of our es- 
timates with the results of the recent applications of the Fade resummation technique to the 
perturbative series for and — ||3^, |3^. It should be stressed that in their analysis 



the authors of Refs. [^, did not consider the light-by-light scattering graphs separately. 
Note also that the coefficients of the corresponding Fade approximants depend from the ln(a;)- 
terms. In spite of the fact that our results for are in qualitative agreement with the results 
of the applications of the Fade resummation method |]35| , [36[| it is interesting to try to under- 
stand the predictive abilities of the Fade resummation methods better. Clearly, this problem 
is connected with the necessity of more detailed understanding of the relations of the Fade 
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results to the ones obtained using the RG-inspired analysis. Note, that the Pade resummation 
methods can face the problem in reproducing the structure of the RG-controlable ln(x)-terms. 
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Appendix 

Using Eq.(14) it is possible to derive the six- loop RG-controUable contributions to the 
diagrams of set (I) : 

d^o^d^^^ ^Be + Ce \n{x) + Dq + Eq + Fq In^(x) + Gq In^(x). (56) 

The expressions for the logarithmic coefficients are 

Ce = 10/30^5 + 8/3iS4 + 6/32^3 + 4/33S2 + 2/34S1 

= 40/3o'S4 + 54/3o/3iS3 + 16/32^2 + 32/30/92^2 + 14/3i/32Si + lW35i 

F, = 80/3o^i32 + ^/3o'/3ii3i 

Ge = 32/3o'Si. (57) 

However, in order to use these expressions in concrete considerations it is necessary to fix 
somehow the value of the five-loop coefficient /34 of the QED /3-function in the OS scheme. 
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